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ABSTRACT 

The prompt emission from GRBs is the brightest electromagnetic emission known yet 
it's origin is not understood. The flux density of individual prompt pulses of a GRB 
can be represented by an analytical expression derived assuming the emission is from 
a thin, ultra-relativistically expanding, uniform, spherical shell over a finite range of 
radii. We present the results of fitting this analytical expression to the lightcurves from 
the four standard Swift BAT energy bands and two standard Swift XRT energy bands 
of 12 bursts. The expression includes the High Latitude Emission (HLE) component 
and the fits provide a rigourous demonstration that the HLE can explain the Rapid 
Decay Phase (RDP) of the prompt emission. The model also accommodates some 
aspects of energy- dependent lag and energy-dependent pulse width, but there are 
features in the data which are not well represented. Some pulses have a hard, narrow 
peak which is not well fitted or a rise and decay which is faster than expected using 
the standard indices derived assuming synchrotron emission from internal shocks, 
although it might be possible to accommodate these features using a different emission 
mechanism within the same overall framework. The luminosity of pulses is correlated 
with the peak energy of the pulse spectrum, Lf oc (E pea k(l + z)) 1 ' 8 , and anti-correlated 
with the time since ejection of the pulse, Lf ex. (T//(l + z))~ 
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INTRODUCTION 



The ability of the Swift s atellite ( GeJrrdset^L 2004h to rapidly and autonomously slew when the Burst Alert Telescope (BAT 
energy range 15-350 keVlBarthelmv et al.l |200Sh ) detects a GRB enables it to point the X-Ray Telescope (XRT energy range 



0.3-10 keV Burrows et al 



( 2005t )) at the target within ~ 100 seconds of the trigger. For many bursts the combination of the 
BAT and XRT provides continuous broadband coverage of the prompt emission through to the onset of t he afterglow. A rapid 



decay phase (RDP) immediately following the prompt emission is observed for the majority of bursts ( Nouseket_aL 20061 ) 



and this decay appears as a smooth continuation of the prompt, both temporally and spectrally ijO'Brien et al.ll2006h . This 



strongly s uggests that t h e dec ay is the tail of the prompt emission and several models have been proposed to account for the 
decay (see IZhang et al. d2007h and references there in) but the most popular is so-called High Latitude Emission (HLE) or 
"naked" GRB emission ( Kumar fc Panaitescu 2000l ) . If the prompt emission from an expanding quasi-spherical shell suddenly 
turns off then radiation from increasingly larger angles relative to the line of sight continue to reach the observer because of 
the added path length (time delay) introduced by the curvature of the shell. The spectrum seen at later (delayed) times is 
subject to a smaller Doppler shift and the flux decays rapidly as the radiation is beamed away from our line of sight, leading 
to a simple predicted relation between the temporal decay index and the spectral index, a — 2 + f3 where F v (t) oc t~ a v~^ ', 
that should hold at late times when t — to ^> At, where to and At are the start time and width of the pulse respectively. 
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Several authors have attempted to demonstrate that the RDP is consistent with the expectations of HLE IjButler fc Kocevski 



2007; IZhang et al.ll2007f ) but simplifying assumptions about to and At have to be made to allow comparison of the model 



predictions and the observations. Usually a parameterized form based on the late time asymptotic behaviour of the HLE 
from a single spike or pulse in the prompt lightcurve is assumed although it is used near the peak of the pulse before the 
asymptotic form is appropriate and the zero time for the powerlaw decay is chosen in an arbitary way. Furthermore, most 
bursts consist of several prompt pul ses which somehow combine to produce the single HLE tail. About 25% of decays have 
no significant spectral evolution and lZhang et al.l (|2007l ) suggest these are dominated by a simple curvature effect while the 
rest show a clear hard-to-soft change as the flux decays which they modelled using an evolving exponential spectrum in which 
the cut-off energy E c = Eo(t/to — l)~ a2 . What is required is a theoretical model that incorporates the pulse profile, the 
evolution of the pulse spectrum, the late HLE tail and summation of the emission from several pulses so that the smooth 
continuation of the flux from the prompt phase in different energy bands can be predicted and compared with the BAT and 
XRT energy-resolved lightcurves. Such a model could be used to fit both prompt pulses seen by the BAT and later soft flares 
seen by the XRT in a consistent way. Our primary aim was to test whether or not the HLE component incorporated into a 
physically realistic model could explain the RDP in a rigourous way but in addition we also confront other aspects of the 
mod el with the data. We h ave used the analytical e xpression fo r the temporal- spec tral profile of i ndivi dual prompt pulses 
from lGenet fc Granotl (j2009r i which is largely based on lSaril (|l998l L [Granotl (|2005l ) and lGranot et al.l {2008). Full details about 
the model i ncluding a compa r ison w ith past attempts to model the RDP with HLE and the properties of the model pulses 
are given in 



Genet fc Granotl (|2009h 



Here we present a first attempt at fitting this model to the Swift BAT and XRT lightcurves of GRBs. We have selected 
f2 objects using the following criteria: early coverage from the XRT including a good measurement of the RDP; relatively 
simple BAT lightcurves; overlapping or near contiguous data from the BAT and XRT giving good coverage of the transition 
from the prompt emission into the RDP; an absence of large X-ray flares which may obscure the RDP. 



2 THE MODEL FUNCTION 



Each individual peak or pulse in the prompt light curve is as sumed to originate from an expanding thin shell ejected from 
the central engine modelled using the emission profile given bv lGenet fc Granotl (|2009h . For the case where the Lorentz factor 



of the expansion, F, is constant with radius (coasting - index m = where F 2 tx R~ m ) the luminosity from within the shell 
in the co moving frame (ind icated by the primed variables) is L' v , = L' (R/ Ro) a S(y' /v' p ). S(x) is assumed to be the Band 
function ( Band et al. 19931 ) and the peak frequency of the vF v spectrum is vL = Vq(R/ Ro) d ■ Under the standard internal 



shock model the index d is linked with the fact that the emission process is synchrotron. Assuming the fast cooling regime 
the peak frequency of the vF v spectrum is v' m oc B' r y^ n where B' is the comoving magnetic field and ■y m is the electron energy 
(Lorentz factor) at the peak of the spectrum. As we are in the coasting phase the strength of shocks is roughly constant, 
so 7 m is constant, and then v' v — v' m tx B' . B' is predominantly normal to the radial direction, so that B' « B /Y oc B for 
m — 0. We expect B oc R~ , hence v' p oc R -1 , and thus d = — 1. More generally, for uniform shells with a roughly constant 
strength of the shocks, both the rate of particles crossing the shock and the average energy per particle are constant with 
radius, implying a constant rate of internal energy generation dE' int /dt' oc R°. For fast cooling this also applies to the total 
comoving luminosity L' ~ vLL' , oc R°, so when L' , evolves with radius so must t/L and we must have d + a = 0. In particular, 
this is true in the case of synchrotron emission where d = — 1 and a = 1. The emission turns on when the shell is at radius Ro 
and turns off at radius Rf — Ro + AR. If m = and d — — 1 we can integrate the comoving luminosity over the equal arrival 
time surface (EATS) giving the number of photons N per unit photon energy E, area A at observed time T as 

' ,A {E ,T^T ej +T ) = P(T-T eh Tf,T rise )B(^^^i] (1) 



dEdAdT" ' - OJ ,- J ,-,. otv _ ^ ^ 

where the pulse profile is 



rp rp \ 1 + 2 / _ \ a + 2 




P(T T e j y Tf,T r ise) — . , ■ . . , 

1 f J \ J-f 

The characteristic times of the pulse Tj and To = Tf — T r ise are the arrival times in the observer frame of the last photon and 
the first photon emitted from the shell, along the line of sight, measured with respect to the the ejection time, T e j. The term 
in square brackets in Equation [2] models the rise in the pulse controlled by the temporal index a and the timescales Trise and 
Tf. The functional form of this rise is the same as in the original formulation but we have included a normalisation such that 
the value is 1 at T — T e j — Tf and have introduced the rise time of the pulse Trise- Figure [T] is a schematic of the pulse profile 
showing the times and timescales. 

B{z) is the Band function with z = (E/E f )((T - T ej )/T f )~ x , such that 
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Figure 1. The pulse profile showing the times T e j and T p j. and the timescales Tf and To = Te — T r . 



B(z) = B nor 



b ^ x e~ z 2O1 -b 2 

b2_1 (6i - b 2 ) bl - b2 e- (bl - b2) z>b 1 -b 2 



(3) 

the characteristic 
= E C (T - T ej ) = 



The parameters of the Band function are the normalisation Bnorm, the low energy spectral index 61, 
energy E c and the high energy spectral index b 2 . The characteristic energy evolves with time, E c (t) 
Ef((T — T e j) /Tf)^ 1 , where Ef is the value of the cut-off energy at time T p k = T e j + Tf. The peak of the vF(y) spectrum 
at this time is E pea k = (61 + l)Ef. T p k marks the start of the RDP and is usually the peak of the pulse as indicated in Fig. 
[1] although if T r i 3e is large, the index d < — 1 and/or the spectrum is hard the pulse m a ximum can occur before T p k. So 
the index d in the present formulation plays a similar role to the index a 2 in lZhang et al.l (|2007l ) although we note that the 
formulation presented here is rigorously valid only for d = — 1, since only in this case does the observed spectrum have the 
same pure Band function form as the comoving emission spectrum. 

The combination of the pulse profile function P(t, Tf,T r i Se ) and the time dependent blue-shift of the spectral profile B(z) 
governed by the index d produce the rise and fall of the pulse. If d = —1 the temporal decay index of the photon count rate 
after t = T — T e j — Tf at frequency v is a v — 2 + /3 V where /3„ is the spectral index of the Band function at that frequency. 
Thus this formulation embodies the well know n closure relation betwe en the temporal and spectral indices from the high 
latitude emission. In the original formulation of iGenet fc Granotl |2009i ) the datum for the characteristic energy of the Band 
function was Eq at time T e j + To. However, we have chosen to use Ef at time T p t because this time marks the maximum of 
the emission when we get the best measurement and estimate of the spectral profile. 

Previous authors have attempted to model the prompt pulses of GRBs, e.g. 



Norris et al.l (|2005l ). The empirical pulse 



profile employed incorporated an exponential rise followed by an exponential decay and the profile was fitted independently 
to each energy band to map any changes of the time constants as a function of energy. The pulse profile used for the current 
work is physically motivated and based on a simple but reasonably realistic model. The rise and fall of the pulse are modelled 
by power laws and the spectral behaviour /evolution is incorporated into the Doppler shift of the spectral profile as a function 
of time. The resulting temporal-spectral profile is fitted simultaneously to all energy bands. 

The prompt emission of a typical burst comprises several pulses and although the spectral-temporal evolution of each 
pulse in the model is governed by the same simple formulation the evolution of a linear superposition of pulses can be complex. 
Each pulse can have a different spectrum (Ef, bi and 62), different time scales (T T iae and Tf and different brighness (B„ or m). 
The current model imposes no restriction on the order or magnitude of the time scales or spectral parameters in the pulse 
sequence and therefore the characteristic energy of the combined spectrum can increase or decrease with time in a complicated 
fashion. 



3 MODEL FITTING 



We have attempted to fit the complete Swift BAT and XRT light curves of the selected sample of 12 GRBs using the above 
m odel for the prompt pul ses. An afterglow component fit to the plateau and final decay of the XRT lightcurve as specified 
bv IWillingale et al.l (|2007l ) was also included so that the transition from the initial decay from the prompt into the plateau 
phase of the afterglow was modelled in a consistent way. The standard energy bands w ere used for th e BAT (15-25 keV, 25-50 
keV, 50-100 keV and 100-350 keV) and XRT (0.3-1.5 keV and 1. 5-10 keV). XSPEC (|Arnaudlll99^ version 12.4.0 was used 
to calculate the expected count rates in these energy bands from the prompt pulse and afterglow components. The number 
of individual count rate predictions required to find a best fit model is very large, ~ 500000 for each pulse or afterglow 
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component. Therefore using XSPEC directly for each evaluation would have been prohibitively time consuming (taking many 
hours to run on a desktop PC) . Instead XSPEC was used to generate a look-up table of count rates over a 3-D grid of spectral 
parameters for the Band function, (61, 61 — 62 and logio(E c )), and the fitting was done by linear interpolation of the predicted 
count rates within the look-up table. A look-up table of 20 x 20 x 20 samples typically takes 20 minutes to generate but the 
fits can then be done in a few seconds. It was found that the Band function model in XSPEC, Grbm, was unable to generate 
spectra with E c < 1.0 keV because there is a hard limit set in the code. Unfortunately late-time spectra from the pulse model 
above can easily reach such low values of characteristic energy so we had to create a local Band function model for XSPEC 
without this restriction. 

The count rates in the low energy XRT band are sensitive to the assumed absorption column density. Fixed values of 
both the Galactic and instrinsic column density were included in the spectral model run in XSPEC using wabs (and zwabs 
if the redshift was known). These column density values were determined separately from spectral fits to the XRT data. For 
the GRBs in the sample the early XRT data are dominated by prompt continuum which may not be a simple power law and 
the intrinsic column density determined using these data might be compromised. We used late time data when available and 
for some bursts with larger intrinsic column density we tried running the fits with different intrinsic column density values to 
see if this influenced the results. The systematic errors introduced by uncertainties in the column densities are smaller than 
the 90% range produced by the fitting so the fits remain essentially the same. 

It was easy to make an initial guess for T p k and T r i se for each pulse. Using this guess a fit was then done letting T p k, 
T r i 3e and Tf vary to find a reasonable fit to the pulse profile and a best fit value for T p k ■ This value of T p k was then fixed in 
all subsequent fitting to find the best fit values of both the spectral and temporal parameters. In some cases the rise of the 
pulse is not well constrained by the data and for these bursts T r i se was fixed. There is in fact some redundancy between the 
T r i 3e and a parameters. A fast rise can be achieved using a = 1 and T r j se small but can also be fitted using a larger T r i 3e 
and a > 1. For most pulses the number and quality of the data points over the rise is such that it is not possible to fit these 
two parameters simultaneously. For fast cooling we expect a + d — and in order to reduce the number of free parameters we 
included this constraint in the fitting. In most cases we fixed a — 1 and hence d = — 1. 

The initial values assumed for Ef and 61 — &2 were 1000 keV and 1 respectively. This specifies a transition to the high 
energy power law at , E = (fei — b2)Ef=10QQ keV. The foi — 62 values derived for bursts observed with BATSE have a considerable 
spread ( Band et al. 1993h but 61 — 62 = 1 is typical. Ef was allowed to float to a lower energy to find the best fit. In a few 
cases the \ 2 was very insensitive to Ef and the best fit energy exceeded 1000 keV. For these pulses Ef was fixed at 1000 
keV. There is some degeneracy /coupling between Ef and the spectral index bi and for most pulses the errors are too large 
to provide independent estimates of these two parameters which determine the spectral shape. In all cases we estimated the 
error range for b\ but we were only able to estimate an error range of Ef for a few pulses. In subsequent analysis we combine 
Ef and b\ to estimate the E pea ^ value and use this single parameter to characterise the hardness of the spectrum. None of the 
fits is able to really constrain 61 — 62 because the BAT response doesn't extend to high enough energies and the total count 
available from individual pulses is usually considerably less than the total fluence. However, 3 soft pulses, pulse 2 in 050724 
and both the pulses in 080805, gave a significantly better fit using a cut-off powerlaw without a high energy tail (61 — &2 — > 00). 
For these pulses we investigated how low 61 — 62 in the Band function could be set before the fit was significantly worse than 
the cut-off powerlaw. In these cases the fixed 61 — 62 used for the fit is effectively an estimate of the lower limit for the change 
in index into the high energy tail of the Band function. 

The fitting was performed using a search for minimum \ 2 an d the resulting minimum \ 2 found along with the number 
of degrees of freedom are listed in Table [T] The product FfT^ +bl for each pulse determines the relative contribution it will 
make to the decaying tail of the prompt emission and all the major pulses visible in the lightcurves were included in the 
modelling. Some faint pulses visible in the data were not included in the model because they were too faint to contribute 
significantly to the RDP. There are features in all the lightcurves which are not well fit by the model and the final \ 2 values 
are not statistically acceptable. Other ways in which the data don't conform to the model can be seen in the plots of the fits 
provided below and will be discussed later. 

Some theoretical/physical constraints were ignored. Pulses were fitted independently with no attempt to introduce cou- 
pling between their parameters. So it is possible that a slow pulse could start before a fast pulse but peak after the same fast 
pulse. In the simple version of the theoretical model with constant P this is clearly not physically possible. We return to this 
point in the discussion below. 



4 FITTING RESULTS 

The Tables [1] and [2] list the fitted parameter values obtained. The sequence of pulses used for each burst are numbered in 
temporal order. When a parameter was allowed to vary the 90% confidence range is given in parentheses with the best fit 
value in the centre. In Table [1] T p k seconds is the time since the trigger of the start of the final decay of the pulse at the 
elapsed time Tf seconds after ejection. This is usually the peak of the pulse but can be a little later. The last column provides 
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the redshift values where available. In Table [2j 61 is the low spectral index of the Band function (i.e. b\ = a + 1 where a 
is the low photon index in the original formulation of the Band function) and 61 — 62 is the difference between the low and 
high spectral indices of the Band function. Ff is the flux in units of 10 -8 ergs cm -2 s _1 over the energy band 0.3 - 350 
keV at time T p k and Ef keV is the characteristic energy of the Band function at that time. The peak energy of the vF(y) 
spectrum at time T p k is given by E pea k = (bi + l)Ef and of course this peak energy evolves with time in the same way as 
the characteristic energy. If bi < —1 then there is no positive peak energy and for these very soft pulses Epeak is probably 
below the XRT energy band. For pulses where the 90% range of 61 encompasses —1 we have not listed an E pea k value. Note, 
we quote the flux over the combined XRT and BAT energy band because all the fits incorporate data from both instruments. 
The last column of Table [5J Lf, is the isotropic luminosity ergs s _1 at time T p k over the bolometic energy band 1-10000 keV. 
This was estimated from the flux Ff using the redshift and assuming a cosmology with Ho = 71 km s _1 Mpc -1 , Q — 0.27 
and A = 0.73 to calculate the luminosity distance cLl- 

We want to stress that the fitted E pea k values listed in Tableware the value of this spectral parameter at a specific time, 
T p k seconds after the BAT trigger. We are able to obtain such values from a Band function fit because we are fitting both 
the spectral profile and the temporal profile of the pulse simultaneously. During a typical pulse the peak energy will range 
from a factor of ~ 1.5, greater than this value, at the start of the rise, to a factor ~ 0.2, smaller than this value, when the tail 
of the pulse disappears into the noise or is obscured by a subsequent pulse. The average pulse spectrum expected from the 
model is a weighted sum or convolution of a series of Band profiles. If the pulse spectrum summed over the main peak were 
fitted with a Band function we would get an average peak energy values, E pea k, which would be similar to but not the same 
as the E pea k values derived by the current technique. Actually, most pulses are too weak to provide sufficient statistics to fit 
a Band function to the average pulse spectrum seen in the Swift BAT. However. iKrimm et al.l (|200E ) were able to fi t spectra 
over relatively small time intervals using a combination of Swift BAT and Suzaku WAM data and lPage et alj (120071 ) provide 
a listing of E pea k values derived for short time intervals of 061121 from cutoff powerlaw fits to Konus-Wmd data. In Table [3] 
we list a small sub-set of the peak energy results which are common between the present work and these studies. The E pea k 
are indeed similar to our E pea k values although they tend to be a little lower because the average pulse spectrum is a little 
softer than that at the peak. When a soft pulse is followed by a hard pulse as is the case for pulses 4 and 5 of 061121 then 
the E pea k value increases as expected. The second time interval for 061121 contains pulse 4 (E pe ak ~ 288) and the start of 
pulse 5 (Epeak = 849) so the E pe ak = 608 is a blend of the two. The third time interval for 061121 contains just pulse 5 and 
is a little harder with E pea k = 621 but not as hard as the peak value for pulse 5. 

The fits for all 12 bursts included are illustrated by the Figures [2] through [13] in date order. The BAT count rates are per 
detector. The observed total count rate is much larger than these values and depends on the position of the burst within the 
BAT field of view. The sub-sections below provide a brief commentary on all the fits including a description of the features 
which are dominating the Ax distribution in each case. We have studied the residuals in detail to see if it likely that the fitted 
parameters and the estimated confidence ranges are compromised by the rather high reduced Chi-squared values. In three 
cases the quoted error range may be too small because of systematic problems with the fitting. These are the 2nd pulse in 
050724 in the 15-25 keV band, the 1st pulse in 060211 A in the 50-100 keV band and the 1st pulse of 080805 in the 15-25 keV 
band. In these pulses the observed flux is systematically higher than the model in the quoted band over a significant fraction 
of the peak. 



4.1 GRB050724 

This is an example of a short burst with an extended emission tail that is just visible in the Swift BAT over a prolonged 
period such that Tgo = 96 seconds. This extended emission is seen as a very broad weak rise in the BAT count rate peaking 
at about 80 seconds and as a very early bright flare in the XRT. This is fitted by the 2nd pulse although it is not possible to 
find a Band function spectrum which fits the BAT and XRT profiles exactly. The fit shown in Fig. [2] has the highest reduced 
Chi-squared value of any of the fits. The model is a little too low over most of the BAT 15-25 keV band and doesn't quite 
match the shape of the tail of the RDP in the XRT 0.3-1.5 keV band. For this essentially short burst the RDP observed is 
entirely due to the extended emission and not the initial short pulse. C oincidentally, this brig ht, soft 2nd pulse gave rise to an 



expanding dust-scattered X-ray halo around this GRB as reported bv lVaughan et al.1 (2006). This burst also has a late flare 



seen in the XRT which is included as the 3rd fitted pulse with T v k = 50000 seconds and Tf = 43000 ± 11000 seconds. 



4.2 GRB050814 

GRB050814 comprises a Fast Rise Exponential Decay (FRED) pulse followed by a weaker 2nd pulse the decay of which was 
tracked well by the XRT. The fitted model can't match the short hard peak of the pulse (seen in the BAT 50-100 keV band) 
which decays rapidly combined with the softer pulse which decays more slowly. The pulse profile is not well fit in the BAT 
15-25 keV band. The rapid decay of the pulse is well fit by the predicted high latitude emission of the model in both XRT 
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Table 1. Temporal parameters for the fitted pulses. T pk , T r i se and Ty are in seconds. A single value indicates a parameter which 
was fixed during fitting. The 90% range and the best fit value a re given for parameters which were allowed to float in the fitting. The 
reds hifts were taken from IProchaska et al] l|2005h . iBloom et al] l|2006l l iFvnbo et al] j2007Mjakobsson et all l|2008l)lThoene et al] j2007l l 
and I Jakobsson et all ifeOOg V 
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Table 2. Spectral parameters for the fitted pulses. Ff 10 — 8 ergs cm~ 2 s _1 over the energy band 0.3-350 keV is the flux at peak, Ef 
keV is the cut-off energy at peak, -E pea fc = (&i + l)i?/ keV is the energy peak of the vF(y) spectrum and Lf ergs s -1 the isotropic 
peak luminosity. A single value indicates a parameter which was fixed during fitting. The 90% range and the best fit value are given for 
parameters which were allowed to float in the fitting. 
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Table 3. Comparison of E pea f s keV at time T p ^ derived in the present work with aver age values, E„ Kn ,k keV, derived for a Band 
function using; combin ed Swift/BAT an d Suzaku Wide-band All-Sky Monitor (WAM) data llKrimm et altoo'sl) and cutoff powerlaw fits 
to Konus- Wind data [jPage et al.ll2007t ). The errors and error ranges quoted are 90% confidence. Both T p j. and Interval are time since 
BAT trigger. 



energy bands except at the transition into the start of the afterglow plateau. This is a good example where the data suggest 
a slower decline before 100 seconds (in the BAT 15-25 keV band) but a faster decline later, in the XRT bands. 



4.3 GRB051001 

A rather weak burst modelled by two pulses. The second pulse is seen in both the BAT and the XRT and the best fit is 
achieved using an Ef = 15 ± 3 keV, which is in between the XRT and BAT band passes. As with other bursts the profile is 
not a good fit in the lowest BAT band 15-25 keV. There is also evidence for curvature of the decay of the 2nd pulse in the 
XRT bands which is not matched by the model. 



4.4 GRB060211A 

This is a rather faint burst seen in the BAT which is reasonably well represented by two pulses except that the pulse profiles 
in the 50-100 keV band have rather fiat tops with a sharp drop at the end which is not emulated by the model. The fit for the 
1st pulse in the 50-100 keV band is particularly poor and it looks as if the spectrum should be harder for this pulse. However, 
if the spectral index is made smaller to accommodate this the decay of the pulse is too shallow, i.e. it is not possible to satisfy 
the closure between the spectral index and temporal index for this pulse. The RDP seen in the XRT is very well matched by 
the model. 



4.5 GRB060814 

GRB060814 was a bright burst consisting of 3 well separated main pulses in the BAT seen with excellent statistics. The best 
fit to the 2nd and 3rd pulse doesn't quite capture the narrow hard peak seen in the BAT 100-350 keV band. The detailed 
shape of the decay of the 2nd pulse is not well fitted and the residuals on this decay contribute significantly to the Chi-squared 
value. We tried varying the index a but this did not improve the fit. We also tried lifting the constraint a + d = and fitting 
a and d independently but this did not improve things either. The 4th pulse is seen in both the BAT and the XRT and it 
proved impossible to fit these data with a single pulse. If we fit to the BAT only then T/ = 28 seconds and the decay falls well 
below the XRT data for t > 150 seconds. Fitting both the BAT and XRT gives T> = 113 seconds but the peaks in the BAT 
are missing and the fit is dominated by the XRT because there are many more data points from the low energy instrument. 
It is not possible to find a spectrum which can produce both the fast, narrow hard pulse and the wider soft decay. 

The fit shown in Fig. [5] consists of a 4th hard pulse which fits the peak seen in the BAT and a slightly later peaking soft 
pulse which fits the decay seen in the XRT. If we calculate the ejection times T e j = T p k — Tf for the 4th and 5th pulses we get 
92 seconds and -102 seconds respectively. Under the simple assumptions of the model this is unphysical since the 5th pulse is 
launched before the 4th pulse (and indeed before the 1st and 2nd pulses!) but the peak is seen after the 4th pulse. 
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4.6 GRB061110A 

This burst comprises a single FRED pulse seen in the BAT followed by a soft pulse seen only in the XRT. The pulse profiles 
and the decay into the afterglow are both well fit by a simple two pulse model. There is a slight mismatch around the transition 
from the rapid decay into the afterglow. 



4.7 GRB061121 

GRB061121 triggered on a well defined precursor which was separated from the main body of the burst by about a minute. 
This gave Swift time to slew onto source and so the XRT was able to captu re several pulses o f the prompt emission alongside 



the BAT. The source exhibited detailed spectral evolution as described bv lPage et al.l (|2007l ). The current model includes 6 
pulses 3 of which were seen in both the BAT and XRT. Because the statistics are good the model is unable to reproduce 
all the features in the data but the overall fit is reasonable. In particular, the cusps of the peaks of pulses 3, 4 and 5 have 
rather large residuals although they are difficult to see on the plot in [8] The fifth and brightest pulse is of particular note. 
The decay from the narrow peak is very rapid but the spectrum is hard, bi = 0.22 ± 0.04 and close to a power law over the 
full BAT and XRT range. It is impossible to fit this decay using the standard model index a — 1.0. The fit shown in Figure [8] 
uses a = 2.4 (and hence d — —2.4). A sixth pulse seen chiefly in the XRT interupts the decay but the overall fit of the high 
latitude emission component from the pulses using this large index value is good. 



4.8 GRB061222A 

A total of 8 pulses was used to model this bright burst and there are still weak features which are not included in the model 
profile. The Chi-squared is dominated by the residuals at the cusps of the pulses 5, 6 and 7 and the transition from the rapid 
decay into the afterglow. 



4.9 GRB070420 

The rather untidy BAT lightcurve of this burst has been fitted using 4 overlapping pulses. The rapid decay seen in the XRT 
1.5-10.0 keV band is not well modelled in the latter stages before it is obscured by the plateau of the afterglow emission at 
~ 200 seconds and this makes a significant contribution to the Chi-squared but the softer decay in the XRT 0.3-1.5 keV band 
is well matched. 



4.10 GRB070621 

This burst is rather different from the others presented here because the lightcurve consists of a series of very similar overlapping 
pulses. The fit shown in Figure [11] comprises 6 pulses. For all these the width T r i se was fixed during the fitting. The BAT data 
are very ragged and clearly more complicated than suggested by the model. However, once again the rapid decay captured 
by the XRT is well fit by the model. The largest residuals occur just before the rise of the 1st pulse. There is clearly some 
emission here which is not included in the model. 



4.11 GRB080229A 

This burst is a good example of a sequence of well defined pulses including a precursor 25 seconds before the main burst. The 
model is unable to fit the hard, narrow peak of pulses 1 (precusor), 2, 3, and 4. The residuals at the cusps of pulses 3 and 4 
a particularly large. 



4.12 GRB080805 

This burst is well represented by just two pulses, the trigger pulse seen just in the BAT and a 2nd soft pulse seen in the XRT. 
The 1st pulse profile is well fit in the BAT 50-100 keV band but the hard peak in the 100-350 keV band is not quite fit and 
there is excess emission seen in the 15-25 keV band which peaks 20 seconds later. We tried many different spectral profiles 
(other than the Band function) to try and match this behaviour without success. The soft tail of this 1st pulse seen at the 
start of the XRT observation in the 1.5-10.0 keV band is also higher than predicted. Both pulses were better fit if Ef was 
allowed to float and the resulting best fit values for this parameter are well constrained. 
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Figure 2. GRB050724 



5 DISCUSSION 



In all cases the fit to the prompt decay from the BAT into the XRT band is reasonably good. In some cases like 050724, 
051001 or 070420 there is a residual curvature in the decay which is not exactly matched by the model but in general the 
summation of the HLE at late times from all the pulses provides a good fit to the RDP. 

There is a generic problem with fitting hard, narrow peaks in conjunction with a soft more slowly decaying tail. The 
fourth pulse in GRB060814 is a prime example. Here it proved impossible to fit the peak in the BAT 50-100 keV band, which 
produces the decay seen in the XRT, using a single pulse. A good fit can be obtained to just the BAT data using a spectral 
index bi = —0.56 but this is far too hard to produce the soft tail seen in the XRT. Alternatively the XRT profile can be fitted 
using bi — —1.09 but this produces no visible flux in the BAT. The fit shown in Figure [6] includes a 5th pulse which fits the 
broad soft tail seen in the XRT. There is a very significant difference between the Tf values of the 4th and 5th pulses in the 
final fit, 38 ± 4 and 252 ± 23 seconds respectively. The required difference in the temporal decay times in the hard BAT and 
soft XRT could not be satisfied using a single spectral profile. A similar sitution exists for the first pulse in GRB080805. The 
profile in the upper three BAT bands is reasonable but this doesn't fit the excess flux seen in the 15-25 keV BAT band or the 
start of the XRT observation in the 1.5-10 keV band before 100 seconds when the second pulse starts to dominate . 

Eoooh 



Despite these problems the model does accommodate an element of energy-dependent lag (jNorris et al 



rowing of the pulses at high energies (jReichart et al.ll2001f ) as illustrated in Fig. 1141 The top panels shows the 1st pulse of 
050814, a typical pulse. The characteristic energy at peak is set to 1000 keV so that the spectrum is very close to a simple 
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Figure 3. GRB050814 



power law over the BAT and XRT energy range at the peak of the pulse. The decay from the peak is a power law in each 
band, slightly steeper for the upper energy bands closer to the characteristic energy. The pulse profile shown in linear scaling 
to the right is very similar in all energy bands. For the 5th pulse in 06f f 21 shown in the middle panels of Fig. [14] the steep 
decay requires a = 2.4 (90% range 2.1-2.8), significantly higher than the standard value of a = 1. Because of curvature in the 
spectrum (the exponential factor in the Band function) the high index also has the effect of shifting the pulse peak to below 
Tf for the harder energy bands. The spectral lag introduced is clearly illustrated in Fig. [14] The best fits for 8 of the pulses 
included Ef as a fitted parameter. For these pulses there is greater spectral curvature because Ej falls within the XRT-BAT 
energy range and consequently the pulse profile changes between the energy bands. The bottom panels of Fig. [14] show an 
example of this; the 1st pulse in 080805. The rise of the pulse is almost identical over all energy bands but the decay is steep 
at high energies and shallow a low energies. The resulting pulse width is a strong function of energy although the peak itself 
remains at the same position. This is not a lag in the strictest sense of the term but it would appear as a lag in any analysis 
involving crosscorrelation since the centroid of the peak clearly shifts to later times as the energy decreases. 

The ability of the model to produce energy dependent lag and/or pulse width is restricted by the requirement for closure 
between the temporal decay index and spectral index. If a pulse is bright in the hard energy band the spectral index tends 
to be close to zero but this in turn tends to produce a slow initial decay of the hard pulse. Conversely, if the pulse is soft 
then the spectral index is large and negative and we expect a fast initial decay in the soft band. It may be possible to fit 
the lag/narrowing seen in pulses like the 1st pulse of 080805 or the hard peak of the 2nd pulse of 060814 using two pulse 
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Figure 4. GRB051001 



components rather than just one. However such component pulses would overlap to a large degree and the parameters of the 
pulses should be coupled in the fitting procedure to prevent the result becoming unphysical. This was not attempted in the 
current work. 

Our model is rather simplistic and doesn't include specific physical details associated with, for example, the internal 
shocks. When two shells collide there will be two emission regions, corresponding to the shocked portions of each shell, and 
they can have the same onset time, T e j + To, the same Ro and T (R ^ Ro) but different emission conditions, E pea k or AR, and 
different peak times, T e j +Tf. We assume the simplest possible history for each shell, that it has coasted at the same V since 
it was ejected, and retains it also when emitting, while in practice multiple collisions are possible and the value of T during 
the emission can be different than before, ft can be increased or decreased due to the collision with the second shell that 
triggers the emission or some deceleration of the front shell due to the sweeping-up of the external medium may also occur. 
The emission may not be caused by internal shocks and may not be synchrotron so the values of indices a = 1 and d = — 1 
might be inappropriate. All these factors will clearly modify the pulse profile and spectrum and the way pulses overlap. 

There are three parameters which were allowed to float in the fitting for all the pulses; the peak flux Ff, the characteristic 
time Tf and the spectral index bi. We have used b\ in conjunction with Ef to calculate the peak energy E pea k which, as 
introduced in lBand et alJ (|1993f ). is a measure of the spectral hardness. The distribution of these parameter values are shown 
in Figure [Pol The top left-hand panel shows Tf vs. the peak time T p k- Somewhat surprisingly there is no strong correlation 
between the position of the peak with respect to the trigger time and length of the pulse although there is a consistent trend 
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Figure 5. GRB060211A 



for the later peaks, T v k > 50 seconds. So within the prompt phase the time scale of the pulses as measured by Tf is not 
dictated by the time since the start of the burst. The sequence of Tf values for 061121 is a good example of this. The top 
right-hand panel shows the correlation between Tf and the peak flux Ff. Here t here is an obvio us tr end which includes the 
late, weak flare seen by the XRT in 050724. Similar correlations were discussed by Lee et al. 1 200oh and Quilligan et al. ( 20021 ) 
although they were considering the pulse width rather than Tf. We demonstra te later in this se c tion t hat T/ is in fact closely 
correlated with pulse width. The first survey of X-ray flares seen in GRBs, IChincarini et al.l (|2007l ) . reports a correlation 
between the flare flux and time of flare since the trigger. Again, this is clearly related to the Ff — Tf correlation although 
we have shown that Tf is not well correlated with the time since trigger for early pulses. The bottom panels show the peak 
energy E pea k, plotted against Tf (left-hand panel) and Ff (right-hand panel). The points for which an error range in Ef was 
calculated are marked in red. For Tf > 10 the peak energy clearly drops with increasing Tf and E pea k is also correlated with 
the peak flux, Ff. 

Table [2] also includes an estimate of the peak luminosity of the pulses for those bursts with a measured redshift. Because 
the peak flux, Ff, is correlated with Tf, as shown in Figure IT5l the luminosity, Lf, is similarly correlated with the characteristic 
pulse time in the source frame, Tf /(1 + z), as shown in Figure lTBI The best fit correlation line plotted is Lf = 1.3 x 10 53 (Tf /(1 + 
ergs s _1 . The most luminous pulse is the 5th pulse in 061121 with T//(l + z) = 0.7 seconds and Lf — 2.6 x 10 53 



2))-2.0±0.2 

ergs s _1 while the least luminous prompt pulses (ignoring the late flare in 050724) were soft and only seen in the XRT, for 
example the 2nd pulse of 050724 with Tf/(1 + z) — 53 seconds and Lf — 3.7 x 10 48 ergs s _1 . A very similar correlation 
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between X-ray flare luminosity and time of flare since the GRB was presented by Lazzati et al. ( 20081 ) . They were considering 
the time since trigger of the peak of the flare, referred to as T p k here, but for late flares seen in the XRT this is very similar 
to Tf. A correlation between Tj and T p k for the later prompt pulses, T pk > 1 00 s, is eviden t in th e top left-hand panel of 
Figure 1151 Assuming that T p k is a good proxy for Tf for the flares analysed by Lazzati et al. ( 20081 ) the values derived here 
extend the relation down to T//(l + z) = 1 second. From this perspective the prompt pulses and late X-ray flares may form a 
continuum consis tent with a common origin associated directly with activity in the central engine. We note that the powerlaw 
index derived bv lLazzati et al. (2008), a = —1.5 ± 0.16, is a little shallower than the value of —2.0 ± 0.19 here although this 
may be due to the small number of pulses/flares included in the analysis. For early prompt pulses there is no correlation 
between T p k and Tf (Figure [15] top left-hand panel) so the extended correlation only makes sense if we consider the time 
since ejection of the pulses not the time since the start of the burst. If this is the case the arguments that the correlation is 
associated with accretion, Lama or alternatively with the spin-down of a magnetar, L oc t~ 2 , are not appropriate. A 
more detailed analysis incorporating a far larger number of prompt pulses and flares is required to see if they really do obey 
a simple powerlaw correlation over the entire range 1 < Tf < 10 4 seconds or whether the prompt pulses and late flares form 
two distinct populations indicating that they arise from different mechanisms. 

A significant correlation of Lf with the peak energy of the spectrum in the source frame, E pea k(l + z), is also shown 
in Figure [TBI The best fit correlation plotted is Lf = 6 .1 x 10 46 (E veak (l + z) keVY - 83±0 ' 16 ergs s" 1 . This correlation is 
similar to the luminosity-peak energy relation reported bv lYonetoku et al.l (j200J) but here we stress that we are considering 
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Figure 7. GRB061110A 



individual pulses and not the average behaviour of bursts. The correlation is between the luminosity at a particular time 
(Tf seconds after ejection when the emission suddenly turns off) and the peak energy at exactly the same time indicating 
a strong physical link between th e these quantities. A spectral-luminosity relation within individual Fermi GRBs has also 
been reported bv lGhirlanda et al. (2009). The correlation they find is consi s tent with the values report ed here and with the 
peak-energy luminosity correlation for long bursts discussed by Nava et al. ( 2008h . Krimm et al. (|2009h provide spectral fits 
with estimates of E pea k for time intervals within individual bursts (see Table [3] and the discussion above) but they concern 
themselves with the E pea k-Ei so relation rather than correlation with the peak luminosity. 

The rise time of the pulse, T r i 3e , was also included in the fitting although in many cases it was fixed because the data 
were not good enough to determine a useful error range. The left panel of Fig. [17] shows a strong correlation between the rise 
time and the characteristic time Tf. As expected T r i 3e is always less than Tf in accordance with the definitions illustrated in 
Fig. [T] This result was not included as a constraint in the fitting but arises from the shape of the pulses fitted. Because of 
this strong correlation the distribution of Trise with respect to the other parameters are very similar to those of Tf shown in 
Fig. [15] The right-hand panel of Fig. 1171 shows the distribution of the ratio T r i Se /Tf with respect to Tf. This ratio is closely 
related to the distance the shell moves during the emission relative to the radius of the shell at the start of the emission 
I Genet fc Granotll2009h . T r i af ,/Tf = AR/(Ro + AR). There is no correlation between this ratio and the other parameters. 

The time from the ejection of the shell to its peak is physically the time when the emission from the shell switches off at 
radius Rf and thus determines the decay time of the peak, which in turn usually dominates the width of the pulse. So while Tf 
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Figure 8. GRB061121 



cannot be measured directly because the ejection of the shell is not observed the width of the pulse largely dictates the value 
of Tf in the fitting and we expect this time to be correlated with the Full Width at Half Maximum (FWHM) of the pulses 
which is directly measureable. The FWHM of the pulses is a complicated function of Tf, T r i Se , the spectral parameters 61, E c 
and bi — and the observed spectral band so we have estimated it by numerical interpolation for the lowest energy band of 
the BAT. The left-hand panel of Fig. H8l shows the FWHM plotted against Tf clearly showing the expected correlation. Note 
that the FWHM is always less or equal to Tf. The right-hand panel shows the FWHM vs. the time of peak since the trigger, 
Tpfe. In this case there is no correlation because the pulse width or Tf are not related to the time since the start of the burst. 
Because of the tight correspondence between pulse width and Tf any correlations involving Tf presented here also hold if 
we substitute the observed FWHM for the fitte d value of Tf . Sim i larly we expect the correlations and analysis involving the 
FWHM of flares given elsewhere, for example in IChincarini et al] ([20071 ). could be interpreted in terms of Tf. 



6 CONCLUSIONS 

To summarise, we have successfully fitted the spectrally resolved light curves of 12 GRBs using an analytical expression for 
pulse profiles derived from the internal shock model incorporating the spectrum and decay indices expected from synchrotron 
emission in the fast cooling regime and a decay tail arising from HLE. A total of 49 pulses describe all the major features 
through the prompt emission into the RDP and, in one case, a late X-ray flare. In all cases the fit to the RDP is good and 
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there is no doubt that the HLE can account for this phase of the emission from GRBs. In many cases the energy dependent 
pulse profile is as predicted by the model showing elements of spectral lag and spectral pulse broadening and indicating that 
the spectral evolution incorporated in the model is in agreement with observations. There are aspects of the data which are 
not in accordance with the simple model. One short, hard pulse in 061121 required index a = 2.4, much higher than the 
value expected for synchrotron emission, a = 1. Several pulses have a hard peak which cannot be adequately fitted. There is 
one pulse, the 4th in 060814, which has a narrow hard peak and broad soft tail which could not be represented by one pulse 
but had to be split into two components. In the present procedure the pulses were fitted independently with no constraints 
between the parameters of individual pulses. Future modelling should include such constraints so that the overlap of the 
pulses remains physically possible. We find that the luminosity of pulses is anti-correlated with the time since ejection for the 
pulses such that bright pulses arise when the time since ejection is small while the dim pulses correspond to long times since 
ejection. The individual pulse luminosity is also correlated with the peak energy of the pulse spectrum consistent with the 
known correlation of peak luminosity for the entire burst with the peak energy of the Tgo spectrum. 
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Figure 14. Pulse profiles as a function of energy band. Top panels: pulse 1 of 050814, a typical pulse. Middle panels: pulse 5 of 061121 
which exhibits energy dependent lag. Bottom panels: pulse 1 of 080805 which exhibits energy dependent pulse width. Left-hand panels: 
the count rate vs. time since T e j with logarithmic scales. The XRT rates have been divided by a factor of 1000 to put them in the same 
range as the BAT rates. Right-hand panels: The count rates have been scaled to a fixed value of 1 at Tf seconds after T e j and axes are 
now linear. In all panels: BAT 100-350 keV solid, BAT 50-100 keV dashed, BAT 25-50 kcV dash-dot, BAT 15-25 keV dotted, XRT 1.5-10 
kcV and 0.3-1.5 kcV both dash-dot-dot-dot. 
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Figure 15. The fitted parameter distributions. Top left-hand panel: The characteristic pulse time, Tt, vs. the position of the peak with 
respect to trigger time, T pk . The late flare in 050724 has been omitted because T pk can be negative and therefore cannot be plotted on 
a logarithmic scale. Top right-hand panel: Tt vs. the peak flux, Ff. Bottom left-hand panel: The peak energy -E pea fe vs. Tt. Bottom 
right-hand panel: The peak energy i? pea fe vs. Ff. The fits for which an error range in Ef was calculated are shown as squares ( red). 
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Figure 16. Left-hand panel: The luminosity at the peak of the pulse Lf ergs s _1 vs. the characteristic pulse time in the source frame, 
Tt/(1 + z). Right-hand panel: The luminosity at the peak of the pulse Lf ergs s -1 vs. the peak energy of the spectrum in the source 
frame, E peak (l + z). 
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Figure 17. The rise time parameter T r i se . Left-hand panel: T, i se vs. Tf. The dotted line indicates equality. Right-hand panel: The ratio 
T r i se /Tf vs. Tf. The points with no error bars for T r i se (i.e. fixed in the fitting) are plotted as squares (red). 




1 10 100 1000 10 4 50 100 150 

T f s T P k s 



Figure 18. The Full Width at Half Maximum of the pulses in the BAT energy band 15-25 keV. Left-hand panel: The FWHM vs. the 
characteristic time TV . The dotted line indicates equality. Right-hand panel: The FWHM vs. the time since trigger of the peak, T p k- 



